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The low reheat temperature at the end of inflation from the gravitino bound constrains the 
creation of heavy Majorana neutrinos associated with models of leptogenesis. However, a detailed 
view of the reheating of the Universe at the end of inflation implies that the maximum temperature 
during reheating, Tmax, can be orders of magnitude higher than the final reheat temperature. This 
then allows for the production of the heavy Majorana neutrinos needed for leptogenesis. We carry 
out the complementary calculation of the gravitino production during reheating and its dependence 
on Tmax. We find that the gravitino abundance generated during reheating for a quartic potential 
is comparable to the standard estimate of the abundance generated after reheating and study its 
consequences for leptogenesis. 
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INTRODUCTION 

It is presumed that the early universe went through a 
period of inflation and then reheated to create the ra- 
diation dominated epoch. If nature is supersymmetric 
then during the process of reheating many supersym- 
metric particles would have been produced, which can 
have important cosmological consequences. In particu- 
lar, the production of gravitinos in the early universe 
and their subsequent evolution including decays has at- 
tracted attention. Stable gravitinos can overdose the uni- 
verse while unstable gravitinos can affect the expansion 
rate of the universe during eras prior to their decay. The 
decay products of unstable gravitinos can also overdose 
the universe or affect light element abundances gener- 
ated during nucleosynthesis. These cosmological conse- 
quences are a function of the gravitino energy density, 
Pq = mQ-riQ, where m^, and tIq are the mass and num- 
ber density of gravitinos. In an inffationary universe, 
Uq is a function of the reheat temperature. Therefore, 
for a fixed mg., often taken to be 0(100 GeV — ITeV), 
cosmological constraints on the energy density of grav- 
itinosprovide an upper bound on the reheat tempera- 
ture % S, i, i, % H, 0, S, i, E El ■ 

The number density of gravitinos is usually obtained 
by considering gravitino production in the radiation dom- 
inated era following reheating, as in Refs. P, 0, [3, 0, H, 
S,0,S,S,[i3,[ii|, and it is presumed that n^, = at the 
beginning of the radiation dominated era. Gravitinos are 
produced through thermal scattering and the final grav- 
itino abundance is found to be proportional to the reheat 
temperature, T^eh- ?reh is the temperature of the ther- 
mal plasma at the beginning of the radiation dominated 
era at trch when the inflaton field has decayed and the 



energy density of the universe is dominated by the infla- 
ton decay products. The cosmological constraints on Uq 
then provide an upper bound on Treh 10^~^ GeV. 

Such an analysis is consistent with the instantaneous 
decay approximation in which one assumes that reheat- 
ing is instantaneous and therefore Ti-oh is the maximum 
temperature during reheating. However, a more detailed 
understanding of (perturbative) reheating indicates that 
during reheating the temperature initially rises to a max- 
imum temperature Tmax and then falls to Tj-eh [13, ■ 
In fact the maximum temperature during the course of 



reheating can be as high as 10^ Tj-eh PJ | 



Earlier works have considered whether sufficient num- 
ber densities of heavy GUT gauge and Higgs bosons, 
or right-handed Majorana neutrinos, required for GUT 
baryogenesis or leptogenesis respectively, can be gener- 
ated with a high Tmax- They find that leptogenesis with 
Majorana neutrinos of mass ~ 10^" GeV is feasible. How- 
ever it is then legitimate to ask if harmful gravitinos are 
also produced during the course of reheating. A priori 
one might expect large production with a Tmax depen- 
dent abundance. This could have serious implications for 
leptogenesis scenarios that invoke large Tmax 13, li , l^l ■ 

In Ref. [l^l we considered this issue for an inffaton 
with a potential of the form V = ^ m^cfP during reheat- 
ing. We find that the abundance of gravitinos produced 
during reheating is 1/3 of that produced in the radia- 
tion dominated epoch. Our estimate for the final grav- 
itino abundance is of the same order as that obtained in 
Ref. [l3| using a numerical analysis (also for a quadratic 
potential). Gravitino production including the reheating 
era contribution has also been obtained numerically in 
Refs. [IMIlSl- In this article we investigate the gravitino 
abundance generated during the course of reheating in 
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inflationary models with a potential V = (A/4)(/)^ during 
reheating. The difference between the two scenarios lies 
in the different equations of state for an oscillating scalar 
field with a (fP and a cj)'^ potential. In the two scenar- 
ios the oscillating scalar field behaves like non-relativistic 
and relativistic matter respectively, i.e., oc R^^, R^^ 
respectively [l^ . This affects the Hubble expansion rate 
in the Boltzmann equation for radiation and gravitinos 
during reheating, and the source term for radiation dur- 
ing reheating, and thus the abundance of gravitinos pro- 
duced by the scattering of the thermalised radiation. 

For the (j)^ potential we find that the gravitino abun- 
dance generated during reheating is about 49% of the 
gravitinos produced in the subsequent radiation domi- 
nated epoch. The gravitino abundance generated during 
reheating is a function of Tmax but, as in the quadratic 
potential case, the abundance can be re-expressed as a 
function of Tj-eh only. Including the contribution from the 
reheating era and then applying the cosmological con- 
straints on the total gravitino abundance lowers the up- 
per bound on T,.oh by a factor of 3/2. This docs not 
have a serious impact on leptogenesis scenarios discussed 
above. 

Our results in Ref. [l^l and in this article are valid 
for chaotic inflation models and for models of inflation 
where one can approximate the inflaton potential during 
reheating by a (fP' or term. However they are not valid 
for a reheating scenario that includes preheating 2l|, |22| . 
23l . |24| | . ^ Gravitino production during prehea ting has 



been considered in Ref. [11, [H, [H, H El 



32, 33, 34, 3 



PERTURBATIVE REHEATING IN 
INFLATIONARY MODELS 

We consider an inflationary model with the inflaton po- 
tential of the form V = (A/4)(/)^ during reheating. ^ The 
inflaton field (p starts oscillating when the inflationary 
epoch ends at a cosmic time t = tosc- While oscillating 
the field cj) decays and the decay products thermalise, ^ 
and thus reheating occurs. ^ Assuming that the universe 



Refs. |25l . |26|| indicate regions of parameter space for a quartic 
potential for which bosonic and fermionic preheating can be sup- 
pressed. Furthermore, certain parameter values for which pre- 
heating is strong produce large non-gaussianities in the CMBR 
and are hence ruled out by WMAP ^27|. 

A quartic term can generate a quadratic term as well. We pre- 
sume that this term, or any other mass term, is small and does 
not dominate till (j> has almost decayed away. 

By decay we refer to dissipation of the energy density of the 
inflaton field due to its coupling with other species. 
We assume that the inflaton products thermalise quickly as dis- 
cussed in Appendix A of Ref. fisll . Refs. [3^ . |37| discuss an 
alternate description of reheating in the context of the MSSM. 
They argue that in the presence of large vevs for flat directions of 
MSSM fields thermalisation slows down. However if the vevs are 



is reheated through the perturbative decay of the inflaton 
field, the reheating picture, in general, can be described 
by O 



Pr 



8ttG 



{Pr + P4>) , 



(1) 

(2) 



where pr and p^ are the energy densities of radiation and 
the inflaton respectively and is the rate of dissipation 
of the inflaton field energy density. Since the equation of 
state for the oscillating inflaton field and for radiation in 
Eq. 12]) is the same we can write the rhs of Eq. ([2]) as 
H'j[Rosc/ RY , where Hj is the Hubble parameter at tosc- 



3 A/pi 



(3) 



where Mj = V^^'^ , Vj being the inflaton energy density 
at tosc- Solving Eq.([2|) then gives 

R = Rose [2Hi{t — tosc) + 1] 



1/2 



(For i> tosc, R^t2.) 

Taking p0 ^ Mf{Rosc/R)^cxp[-T4,{t - to 
^ we then get 

3 „2 „2 1 - e-r*(t-to=.) 

Pr — -'iJpi-Wt n 

StT ' [2H,{t-tosc) + l]' 



(4) 
in Eq. 

(5) 



However to simplify our subsequent analysis, we ignore 
the change in p^ due to decay in Eq. ([T]), which is valid 



till t < i,eh 



solution of Eq. ([T]) is given by 



Then 



Mf{Rosc/Rr and the 



^ _ 3 2 p 7:t2 t ^osc 



R 



327r ' ^ '^'Vi?osc 
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1 - 



R 

^OSC 



(6) 



(7) 



From Eq. ([7| we see that during reheating the energy 
density initially increases to a maximum value 



4 V 327r ^ ^ 



(8) 



at i?max — V^R 

reheating is then 
T — 

max — 



The maximum temperature during 



0.6g."'/^F^Mpi)^''*Af 



1/4 ,,fl/2 



(9) 



small (< 10~^Mpi) then there is no effect on reheating. Further- 
more the condensates of MSSM fields can fragment into Q-balls 
and the vevs can vanish in large parts of the universe [38| . These 
are the cases we would be considering. 

We will assume for now that our analysis below is valid till tj-^^^ 
and will later discuss this assumption. 
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Subsequently the temperature falls as (for t :$> 

imax) until the inflaton decays at treh- Once the final 
decay products of cj) thermalise with each other through 
sufficient interactions the radiation density becomes 



30 



(10) 



From equations ([6]) and (|10p . and assuming tj-^h ^ 
we get the reheating temperature 



rch 



(11) 



In the following we examine the production of graviti- 
nos during reheating, i.e., from iosc to t^ch, and during 
the subsequent radiation dominated era after treh, and 
discuss its consequences. 



with the gauge groups. ci,2,3 are 11, 27 and 72 and ki 2,3 
are 1.266, 1.312, 1.271 respectively (Table 1 of Ref. 0]'). 
The above expression includes corrections to earlier ex- 
pressions for the cross section for gravitino production 
in Refs. (i^l and 17 1. Using the one loop /3- function of 
MSSM, the solution of the renormalization group equa- 
tion for the gauge coupling constants is given by 



1-1/2 



(15) 



with 61 = 11, 62 = 1, 63 = —3. To obtain a conservative 
estimate of the gravitino abundance we take Mi —s- as 
in Ref. [3- 



Gravitino production during reheating 



GRAVITINO PRODUCTION 

Gravitinos arc produced by the scattering of the in- 
flaton decay products; a list of processes is provided in, 
for example. Tables 1 in Refs. [l3|- The Boltzmann 
equation for gravitinos is given by 



(12) 



where n = {({3)/7t'^)T^ is the number density of scat- 
terers (C(3) = 1.20206.. is the Riemann zcta function 
of 3), Stot is the total scattering cross section for grav- 
itino production, v is the relative velocity of the incoming 
particles, and (...) refers to thermal averaging. Since the 
gravitino lifetime is 10^~^(100 GeV/m(5)s [5| decays are 
not relevant during the epoch of gravitino production for 
gravitinos of mass 10^~^ GeV. Hence we have not in- 
cluded the gravitino decay term in Eq. I|12p . We may 
now re-express this equation as 



+ 3Hn^ = (S, 



\v\)r 



(13) 



(keeping in mind that T passes through zero at Tj^ax)- 
The cross section (Stotbl) is given by j39| 



(Stot|w[ 



a 

M2 
1 



3tt 



3 



16C(3) ^ 



3m| 



X a g, 



2 In 



(14) 



where i = 1,2,3 refers to the three gauge groups 
U{1)y,SU{2)l and SU{3)c respectively. M = 
Mpi/\/87r ~ 2.4 X 10^® GeV is the reduced Planck mass. 
Mi are the gaugino masses. gi{T) are the gauge cou- 
pling constants, while Ci and ki are constants associated 



Since T is zero at Tmax we solve Eq. (fT5)) from ^osc to 
imax and from tmax to treh Separately. In order to solve 
Eq. ([Ill) we need f and H as functions of T. Eq. ^ 
implies that 



R'^T^ — AR^ R'^ + AR 



4 

osc 







(16) 



where with some algebra one can show that A = AT^^ 
Then in the epoch tosc l£ t < tmax we have 



R'^/R; 



OSC 



2 

osc 



i-(i-4rVA)2 

2T^/A 



1 + {1-AT^/A)l 



2T^/A 



(17) 



(18) 



Eq. (US]) imphes that for t > i,nax, T ^ R~^ (not ~ R^^ 
as inflaton decay is a source of radiation) . We now define 
a dimensionless variable 



4^4 j,A 

a; = 1 T- = 1 r- 



(19) 



Then the Hubble expansion parameter H = iJ/(i?osc/ R)^ 
can be rewritten as a function of T or x as 

H^^Hi{l + x) for W<i<tmax (20) 

and 



H=-Hi{\-x) for W<i<iroh. (21) 

For T we differentiate p,.(T) = {■k'^ /3Q)g^T* and p,.(i?) 
in Eq. ([7]) with respect to time, equate the results and 

get 

TT^g* 4T^ dR Apr dR ^ ' 
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Now using Eqs. (|17ll8p . (fT9l) and (|20l2ip the above 
equation can be recast as 



T = 



rpA TT 

2T3 



.t(1 + for iosc < t < i max (23) 



and 



T 



max I 

2r3 



(-.t)(1-x)2 for t,„ax < i < trch . (24) 



Wc now solve the Boltzmann equation Eq. ([T5)) in the 

two regimes iosc < ^ < ^max and tmax <t< trch- 



Epoch: tosc <t<tr, 



Using Eqs. ([20]) and (|23| Eq. ([TS]) can be written as 



driQ 



dx (1 + x) 



-G = rf2^V^y^ (25) 



(1+X)2 



where 



Now we define y = 1+x, so that Eq. ([^5]) can be rewritten 
as 

(2-y)3/2 



1/2 



(27) 



Solving Eq. (|?7|) from ?/osc = 2 to y and assuming that 
n-ciVosc) = we get the gravitino abundance as [41| 



3 



2(2 -y)- 6 



x/2 



In 



(28) 



Thus at y = ?/max = 1, which corresponds to T = Tmax, 
we get the gravitino abundance to be 



'^^(ymax) = d2 



-4 ^ In 

V2 



1 - \/2 



l + \/2 



-0.26 da 



(29) 



Solving Eq. (13 1| from ymax to y we get the gravitino 
abundance as |4l| 



nciv) = '^G(ymax) 



y 



ymax 



+ d,y^'^ 



4 + ^ In 
%/2 



3/2 



1- V2 



V2 



1 + V2 
2(l + i)v^ 

^/2~i-^/2 



V2~y + \/2 



(32) 



Now using Eq. (j29p and letting j/max = 1 and ?/ = j/roh 
the above equation can be written as 



«G(yrch) = cf2 2;reh 



4^- 



-2(1 + ) v/2 - yrch 



V2 - 2/rch - \/2 



V2 - 2/rch + \/2 



(33) 



Using J/rch ~ (l/2)(r,.ch/Tniax)^ for Tioh < r,nax WC gCt 

the gravitino abundance at Trch as 



"0(71011) 



23/2 



(^rch/^max) 



_4^/2 ('^ 

J rch 

8T4 



3 Tt 



8 



1 7^8 

^ -'rcl 

16 



V2 '\l6T4,,-r4;^ 



(34) 



For Ti-eh ^ Tmax the dominant contribution comes from 
the (T'lnax/T'rch)'' term in the square brackets. Then we 
can approximate the gravitational abundance at Trch as 



"G(T'roh) ^ ^575- (Trch /Tir 



-4V2 



Hr 



(35) 



Gravitino production in the radiation dominated era 



Epoch: tmax < t < trch 



Using Eqs. and ([M]) Eq. can be written as 



dn^ di 



dx (1 — a;) 



(l-x^)3/2 



dng di , (2 - zj)3/2 

dy y ^ y^i^ 



(30) 



where di and d2 are given by Eq. p6|) . Defining y = 1 — a; 
Eq. can be rewritten as 



(31) 



After the inflaton field decays at i^ch the universe en- 
ters the radiation dominated era. Unlike the reheating 
era during which the entropy continuously increases, in 
the radiation dominated era the total entropy remain 
constant (except for epochs of out-of-equilibrium decays). 
Therefore it is useful to express the abundance of any 
species « as = ni/s. where rii is the number density of 
the species i in a physical volume and s is the entropy 
density given by 



(36) 
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FIG. 1: The production of gravitinos is shown from fosc to 
trch through fmax as a function of temperature for values men- 
tioned in the Discussion. The temperature rises from at tosc 
to Tmax at tmax and then falls to Trch at trch- The (nor- 
malised) gravitino number density Ri = nQ(r)/ng(Tmax) 
rises from to its maximum value at tMAx > tmax, when 
Tmax — O.STmax, and then decreases till troh [red solid line]. 
The (normalised) gravitino number density per comoving vol- 
ume R2 = ficiT) / nciTinax) is also plotted [blue dashed line]. 
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FIG. 2: Yq — Uq/s generated during the reheating era, the 
radiation dominated era, and the sum of contributions from 
both eras are shown as a function of the temperature T for 
t > tmax. Treh and Ml are chosen to be 10^ GeV and 10^*^ GeV 
respectively, and so Tmax ~ 1.3 x 10^^ GeV. Since Yq in both 
eras is largely generated close to Treh, a is evaluated at T^ch- 
The final value of V^. is ~ 3 x 10"". 



We take 5* = 228.75 in the MSSM for the temperature 
range of interest. One can now re-express Eq. ((T3|) as 



• f]Y- 



(37) 



To obtain T we use the temperature-time relation for the 
radiation dominated era, namely, 

T = Trch r : (38) 



[2i/reh(t-troh) + l]- 



where 



90 Mpi 

(For t > trch, T ^ t^i.) Therefore T is given by 



90 / M 



(39) 



(40) 



Then 



dT 



90 [ 45 \ /^x /C(3) 



(41) 

Assuming a to be independent of temperature and in- 
tegrating the above equation from Tj-ch to Tf, the final 
temperature, we get the total gravitino abundance at Tf 
to be 



Y^{Ti)^Y^{T,,^,)+Yl-^iTf) 



(42) 



where 



1/2 



/ 45 



is the gravitino abundance produced in the radiation 

dominated era. We have used g*reh in the expression 

for Y'^'^'^ and ignored the variation of 5* with tempera- 

^ . . . . . 
ture. This is justified since most of the gravitinos are 

generated close to T^oh- Using Eqs. ((35)) and ([36l) 



^G(rreh) = T72 



2a /C(3) 



A/2 V TT^ 



45 



2TT^g^, 



rp4 



(44) 



This contribution is usually neglected while obtaining the 
gravitino bound. As we see below this is comparable with 
the second term in Eq. (gl]). Using Eqs. (gS]) and (gH) 
in Eq. gH) we get 



A/ 



45 



.g*roli7r" 



27r2g, 
1/2 



rp4 



(45) 



where we have used Tf ^ Tich- Relating Tmax to Trch 
from Eqs. ([3]), ^ and pTjl the total gravitino abundance 
is then given by 



3aT^ch f C(3) 



45 



O 2 3/2 



(46) 



A/ V 
[0.49 + 1.0] , 

where we have used _g*rch in the expressions for T„ 
DISCUSSION 



The detailed dynamics of gravitino production from 
tosc to troh is shown in Fig. ([T|) as a function of the 
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temperature using Eqs. p8|) and ([32)) . We normalise 
the gravitino number density with respect to the value 
at Tmax- ?rch and Mi are chosen to be 10^ GcV and 
10^^ GeV respectively, and so T^ax ~ 1-3 x lO^^ Qgy. a is 
treated as constant and evaluated at Tioh since most grav- 
itinos are produced near Trch- a is 15.1, using gi{Mz) ob- 
tained from aEAiiMz) = 1/128, sin^ 0w(Mz) = 0.231, 
asiMz) = 0.119, and Mz = 91.2 GeV 

It can be seen in Fig. ([1]) that during reheating 
the number density of gravitinos monotonically increases 
from tosc to a time ^max > imax, where iiviAX corresponds 
to a temperature Tmax — O.STmax and nQ{tMAx) ~ 
3n(3(imax)- Subsequently it decreases till troh to a value 
7.7 (Trch/rmax)^ »^G(imax)- Taking the number density 
per comoving volume as nc{T) = {tIqR^) / R^^^ we also 
plot no (normalised to the value at Tmax) by using Eqs. 
([TTll. mi), jMl) and ([32]). From the plot ofnaiT) in Fig. 
^ it can be seen that most gravitinos during reheating 
are produced close to T^eh- 

In Fig. ([2]) we show the contribution to the gravitino 
abundance from the reheating era and from the subse- 
quent radiation dominated era, and the sum of these con- 
tributions. From Eq. (^5)) it is clear that the gravitino 
production during the reheating era is almost half of that 
during the radiation dominated era even though a priori 
one would not have expected the gravitino production in 
both these eras to be similar. While the gravitino abun- 
dance generated during reheating is a function of T„iax it 
is interesting that it can be re-cxprcssed as independent 
of T,„ax, and as a function of only Tioh. Moreover, the 
contribution to Y^^ from the reheating era is linearly pro- 
portional to Treh, as it is for the radiation dominated era. 
These results are similar to those obtained in Ref. [3]. 
The linear dependence on Tj-eh makes it simple to revise 
the constraints on Treh based on the upper limit on the 
gravitino abundance - the upper bound on Trch is low- 
ered by a factor of 3/2. Since Tmax oc VTrcin 7i„ax is 
not affected much. Therefore models of leptogenesis that 
invoke a large T,„ax to create heavy Majorana neutrinos 
are not significantly impacted. 

Above we partly ignored inflaton decay in our analysis, 
i.e., we did not include the effect of cxp[~T^{t — tosc)] 
in in Eq. ([1]). One might be concerned that this 
will lead to inaccuracies close to irch when most of 
the gravitinos arc produced. However if one writes 

/90 ~ i?~'^exp(-r0t) t~'^ C-K.p{-T^t) for t > tmax 

then P4,/ P(f, — — 2/t — F^. Therefore even till close to 
trch ~ p^ decreases primarily due to the expansion 
of the universe. Furthermore, if we follow the value of 
Yq = tIq/s we find (from Eq. (|44|) ) that near t^^h it in- 
creases as T-^ ^ ^ n. AitKi O.IF^^ 56% of our 
estimate of Yg, is already generated while decay has led to 
a reduction in p^ of only 9%. Keeping in mind the above 
arguments, we expect that the error in our estimate of 
Yq will not be large. A more accurate estimate will re- 
quire a numerical analysis. (For the quadratic potential. 



Ref. [ll] (analytic) and Ref. [l3| (numerical) obtained a 
gravitino abundance of 1.9 x 10"^'^ and 1.5 x 10"^'^ re- 
spectively for Treh, as defined in Ref. set to 10^ GeV. 



CONCLUSION 

In conclusion, in this article we have calculated the 
gravitino abundance generated during reheating for an 
inflationary model with a quartic potential during re- 
heating. We find that the gravitino abundance generated 
during reheating is a function of the largest tempera- 
ture during reheating. However it can be re-expressed in 
terms of the reheat temperature only and we find that it 
is linearly proportional to the the reheat temperature, as 
in the standard calculation of gravitinos produced in the 
radiation dominated era after reheating. Furthermore, 
we find that this abundance is 49% of the abundance of 
gravitinos generated in the radiation dominated era. This 
lowers the upper bound on the reheat temperature by a 
factor of 3/2. However this does not significant change 
the viability of leptogenesis scenarios. 
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